MIFP-09-19, NSF-KITP-09-52, arXiv: 0904.3773 



Mirage in the Sky: Non-thermal Dark Matter, Gravitino Problem, and Cosmic Ray 

Anomalies 

Bhaskar Dutta^Q Louis Leblond^'^ and Kuver Sinha^ 
^ George P. & Cynthia W. Mitchell Institute for Fundamental Physics, 
Texas A&M University, College Station, TX 77843-4242 
^ Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA, 93106, USA 

Recent anomalies in cosmic rays could be due to dark matter annihilation in our galaxy. In 
order to get the required large cross-section to explain the data while still obtaining the right relic 
density, we rely on a non standard thermal history between dark matter freeze-out and Big-Bang 
Nucleosynthesis (BEN). We show that through a reheating phase from the decay of a heavy moduli or 
even the gravitino, we can produce the right relic density of dark matter if its self-annihilation cross- 
section is large enough. In addition to fitting the recent data, this scenario solves the cosmological 
moduli and gravitino problems. We illustrate this mechanism with a specific example in the context 
of U(\)b-l extended MSSM where supersymmetry is broken via mirage mediation. These string 
motivated models naturally contain heavy moduli decaying to the gravitino, whose subsequent decay 
to the LSP can reheat the universe at a low temperature. The right-handed sneutrino and the B — L 
gaugino can both be viable dark matter candidates with large cross-section. They are leptophilic 
because of B — L charges. We also show that it is possible to distinguish the non-thermal from 
the thermal scenario (using Sommerfeld enhancement) in direct detection experiments for certain 
regions of parameter space. 



I. INTRODUCTION 

Recent data from PAMELA shows an excess of 
positrons at energies in the 10-100 GeV range W, with 
data expected up to ~ 270 GeV. No excess of anti- 
proton flux is observed [2] . There is also new data from 
ATIC where one observes excess in e+ -I- e~ spectrum 
with a peak around 600 GeV [3]. PPB-BETS g], also 
reports excess in the -I- e~ energy spectrum between 
500-800 GeV. While there could be astrophysical expla- 
nations for these anomalies [5], it is also possible that 
these are among the first signals of dark matter annihi- 
lation. 

If PAMELA is explained by a WIMP dark matter, the 
data leads us to the following three broad characteristics 
for this particle: It must be heavier than ~ 100 GeV, 
it must be leptophilic and it must have a large cross- 
section today [B] . The first property is needed to explain 
the high energy positron detected while we need to have 
final states of dark matter annihilation predominantly 
to be leptons in order to not overproduce anti-protons 
[Zl [8]. Both of these properties can be easily arranged 
in a model dependent way. The necessarily large cross- 
section on the other hand is harder to fiddle with since it 
is directly constrained by the relic density. The thermally 
produced relic density is given by 

ncoM = 0.23 1^) j ■ W 

An interesting proposal is to enhance the cross-section at 
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low-velocity through the Sommerfeld effect [51 whereas 
a light boson provides an attractive potential that en- 
hances the cross-section when the dark matter is non- 
relativistic [71 [TU] (also see [TTJ [T^] for explanations using 
Sommerfield effect). In order to generate the right en- 
hancement factor one needs to fix the ratio of the dark 
matter mass and the new light boson mass (or different 
parameters of the model) to a high degree of accuracy. 

In this paper, we will be interested in a second alterna- 
tive where we have a non-standard thermal history and 
Eq. ([ij is modified. This possibility has already been ex- 
plored as a possible explanation for PAMELA in [T3j [14] 
(see also |inj)- Such a non standard thermal history is 
very well-motivated from physics beyond the standard 
model and, of course, we have no direct evidence that the 
universe is radiation dominated at temperature above the 
Big-Bang Nucleosynthesis (BBN) temperature. If there 
is a phase of matter or dark energy domination prior to 
BBN, the only way to connect to radiation (whose en- 
ergy density decays faster than either) is through a re- 
heating process [TB]. If the reheating temperature (T^) 
is in between the freeze-out temperature (Tf) and the 
BBN temperature, we will respect all current astrophys- 
ical constraints while the entropy produced at reheating 
will naturally dilute the relic density produced at freeze- 
out. If dark matter is produced non-thermally at the 
time of reheating, larger annihilation cross-sections are 
needed to obtain the right relic density. Similarly, we 
could have a phase dominated by a fiuid whose energy 
density decays faster than radiation (e.g. kination [T7'). 
In this case, there is no reheating phase but the Hub- 
ble expansion is faster than usual and the dark matter 
candidate must have a stronger coupling in order to get 
the right relic density. We therefore see that almost any 
thermal history other than radiation domination would 
require a dark matter candidate with larger annihilation 
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cross-section to obtain the correct relic density. 

To illustrate this non-standard thermal history, we 
work with the specific example of a phase of matter dom- 
ination before BBN. We consider the case where the mat- 
ter component is a scalar field coherently oscillating (a 
cosmological modulus) and also the case where the mat- 
ter component is the gravitino. The former case is well 
motivated from string theory where there are many flat 
directions (moduli) that acquire masses from supersym- 
metry breaking. 

A second purpose of this paper is to explore the im- 
plications of such cosmological enhancement of the dark 
matter annihilation cross section on the cosmological 
moduli and gravitino problems. Generic string mod- 
uli of masses around the electroweak scale decay (and 
reheat the universe) after BBN, ruining its successful 
predictions. A standard solution to this problem is to 
take the moduli to be heavy (at least 20 TeV), thus en- 
abling them to decay before BBN. The decay of such 
moduli primarily produces gauge bosons, gauginos and 
dark matter; however, there is also gravitini production 
which is generically unsuppressed, with a branching ratio 
of around 0.01 [T5] (this can be avoided in special set-ups 
[19p. The production of gravitinos is again problematic 
and they need to be heavy enough to decay before BBN. 
Moreover, the gravitinos decay to the LSP, and there are 
strong bounds to avoid overproduction of dark matter. 
Avoiding overproduction can cause the gravitino mass to 
be high (around 1000 TeV) if the annihilation cross sec- 
tion of LSP is at the canonical value. This has disastrous 
implications for low-energy supersymmetric model build- 
ing. Superconformal anomaly mediated contributions to 
the soft masses push the low energy superparticle spec- 
trum into the 100 TeV region. 

We will show that a large enhancement of the cross- 
section (~ 10^) and heavier dark matter (TeV scale) can 
naturally ease the bounds on gravitino mass coming from 
overproduction of LSP. The gravitino becomes a legiti- 
mate candidate (together with cosmological moduli) to 
act as the decaying particle that non-thermally produces 
the right relic density of LSP. We thus solve the grav- 
itino problem by having it to decay prior to BBN, while 
in order not to overproduce the LSP we need the large 
cross-section required by PAMELA! 

Non-thermal dark matter production and the moduli 
and gravitino problems have usually been studied in the 
context of a Wino LSP, which arises in models of anomaly 
mediation |20) , simple realizations of split supersymme- 
try, and in the context of the G2-MSSM |5T1. 

In view of recent data, the high energy positron excess 
reported by PAMELA is difficult to fit with a LSP in the 
Wino mass range, unless non-standard assumptions are 
made for the distribution of dark matter and the propa- 
gation of cosmic rays (this has been studied extensively in 
[Uj). Moreover, while it is possible that the anti-proton 
data suffers from theoretical uncertainties in cosmic ray 
propagation, taken at face value such data appears not 
to prefer a Wino LSP. 



The point of view we will take in this paper is that the 
above cosmological scenario can work in a U{1)b~l ex- 
tension of the MSSM, in the setting of mirage mediation. 
From a model-building perspective, the fact that a non- 
minimal model eases bounds is perhaps not entirely sur- 
prising. However, as we will show, this particular exten- 
sion (already well-motivated by non-zero neutrino mass) 
has a TeV-scale leptophilic LSP (the right-handed sneu- 
trino or the B — L gaugino) with large cross section. The 
string inspired models of mirage mediation also solve the 
tachyonic slepton problem of anomaly mediation. 

In comparison to the usual thermal production of dark 
matter, an enhancement factor given by the ratio of the 
freeze-out temperature to the reheat temperature gets 
generated in this scenario. After solving the moduli and 
the gravition problem, we will show that the enhance- 
ment factor is of order ^ ~ 10^~* which is in the right 
range to explain the cosmic ray puzzle. 

We note that for mirage mediation in the MSSM, the 
LSP is primarily the Bino, which is unacceptable in light 
of the PAMELA data while our dark matter candidates 
in the U{1)b-l extension can fare better. Thus, apart 
from solving the moduli/gravitino problem, our model 
connects (possible) indirect observation of dark matter 
with string inspired phenomenology. We also show that 
in the case of a right-handed sneutrino, it is possible to 
distinguish the non-thermal from the thermal scenarios 
in direct detection experiments. 

In section §2 we give details of the cosmological en- 
hancement, including the gravitino problem. In section 
§3, we work out the example of the mirage mediated B—L 
extension, which provides a concrete model where the 
above cosmological history can occur. We conclude in 
section §4. 

II. COSMOLOGICAL ENHANCEMENT 

Here we will work out the cosmological enhancement 
of the cross section caused by the reheating of a heavy 
modulus before BBN (the story is similar for the grav- 
itino decay). The key phenomenon here is the low tem- 
perature of reheat (LTR) [16 which could also come from 
a phase of dark energy dominance such as low scale infla- 
tion or thermal inflation instead of a phase dominated by 
matter. Cosmological moduli can come to dominate the 
energy density of the universe if they are displaced from 
the minimum of their potential. The equation of motion 
of a scalar field with gravitational strength decay rate in 
a FRW background is 

'(i)+{'iH + V4,)4> + V' = Q . (2) 

After infiation, the initial vev of the field {(t>in) is dis- 
placed from its zero temperature minimum by some 
amount M (say inflationary scale). At early time, H > 
m^, the friction term dominates over the potential and 
the field is frozen at its initial value (f> = (pin- The uni- 
verse is then radiation dominated until t ^ 'TiT^ (at a 
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temperature Tin ~ y/iT^pMp) at which point the field 
will start oscillating around its minimum. These coher- 
ent oscillations of a scalar field lead to large occupation 
numbers and the end result is a Bose-Einstein conden- 
sate which behaves like matter. The initial energy den- 
sity (TO00f„) will increase compared to radiation and it 
will eventually dominate until the modulus decays and 
reheats. If the modulus decays after BBN, the energy re- 
leased will photo-dissociate the newly built nucleus [H] 
which is the crux of the cosmological moduli problem. In 
the following, we will take the cross-section of the mod- 
ulus to ordinary matter and to dark matter to be Planck 
suppressed 



2t: M„ 



(3) 



where we take c ~ 1 and Mp — 2.4 x lO^'^GeV is the 
reduced Planck mass. In the approximation of sudden 
decay the reheating temperature can be defined by tak- 
ing the lifetime of the modulus (P^^^) to be equal to the 
expansion rate at the time of reheating t = Right 
after reheating (at T = Tj.) the universe is radiation dom- 



inated with H : 



90 A/„ 



where g^, counts the number 



of degrees of freedom^ . Since our temperature of reheat 
will be always be around 10 — 100 MeV, has the usual 
value of = 10.75. The temperature of reheat is then 



1/2 3/2 



T 



.1/2 



10.75 

9* 



100 TeV 



3/2 



6.37 MeV. (4) 



We can compute the relic density of dark matter using 
the Boltzmann equations for the modulus (j), the dark 
matter candidate X and radiation R [20j 



dt 
dpR 
dt 



dn 



X 



dt 



-3Hp^ - P^ 
-AHpii + (to0 - NLSP'mx)T^n^ 
-{(7v)2mx 

-3Hnx + NLSpT^cpnx 



-{av) n\ - (n^f 



(5) 



(6) 



(7) 



where N^sp is the average number of LSP particles 
produced by the decay of one modulus and we have 
taken the energy of both (p and X to be simply given 
by their masses, i.e. they are non-relativistic. Differ- 
ent initial conditions at freeze-out are possible and we 



^ This relation is modified by 0(1) factor because of the non- 
standard cosmology |16l with the massive modulus. 



will look at the case where the universe is dominated 
by the moduli prior to freeze-out^. For high masses and 
strong enough interactions, the dark matter candidate 
will be non-relativistic at the time of freeze-out (with 



3/2 



) and it will have time to 
reach equilibrium before reheating occurs. Dark mat- 
ter freeze-out occur when the annihilation rate is equal 
to the rate of expansion 



P 



X 



n'^^iTf) (av) = HiTf) 



(8) 



up to the fact that the thermal history is changed due to 
the presence of the decaying moduli. This change is rel- 
atively minor and the new freeze-out temperature Tj?'^'" 

is very close to the usual one at Tj''' ^ mx/20 which is 
what we will be using. After freeze-out, reheating occurs 
and the entropy production will dilute the initial density 
of dark matter by a total factor of /Tj which can be as 
much as 10^^^. To first order, we can therefore neglect 
the initial density of dark matter coming from freeze-out 
and instead just focus on the component produced non- 
thermally from the decay of the heavy moduli. 

To compute the non-thermally produced dark matter, 
we can use the fact that there is an attractor solution to 
the Boltzmann equations. The idea is that if dark matter 
is overproduced by the moduli (compared to what one 
has for a usual freeze-out) they will quickly annihilate 
back into radiation. Therefore the maximal density of 
dark matter is given by the same condition we used before 
(Eq. [8| but now at a lower temperature 

n'^j^iTr) (av) = H{Tr) . (9) 

The non-thermal density of dark matter scales like 



nx{Tr) H{Tr) 1 



s{Tr) 



(10) 



where s = "^j^g^T^ is the entropy density. This should 
be compared to the usual thermal freeze-out density 



nl{Tr) 



1 



(11) 



Hence the non-thermal production is enhanced compared 
to the usual thermal one by a factor 



NT 
X 



T 



X 



T 

J. r 



0.23 I 

(av) J Tr 



and we must increase the cross-section accordingly to ex- 
plain the data (by the factor Tf/Tr). 



^ Demanding that the universe is matter dominated prior to freeze- 
out imposes the following condition on the initial value of 4>, 

<f>in > fl/^ ^^1/4 which is around lo'-'-^^-'-^) GcV for the numbers 

considered in this paper. This is well below the expected value 

of ~ Mp. 
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In the case, where the modulus has a small branching 
ratio to the dark matter particle (small Nlsp) or in the 
case where the modulus does not dominate the energy 
density before it decays, one may not reach the attractor 
solution. In this case, one can show that the dark matter 
abundance (or yield) Yx{T) = is just given by the 
abundance of the modulus at reheating times the branch- 
ing ratio Brjj^xYrj,. If one work with number density in- 
stead of abundance, we should use the average number 
of particle produced N^sp instead of the branching ratio 
as in in]. 

Therefore the abundance of LSP is the minimum 

(13) 

The first possibility represents the case where not enough 
dark matter is produced for self annihilation to start 
while the second is the attractor solution described 
above. At this point the mass of the modulus is a free 
parameter and we can tune it to get any temperature of 
reheat desired while tuning the cross-section accordingly 
to get the right relic density . For a very heavy modulus, 
an important worry is that gravitino will be produced in 
the reheating phase creating a new problem (or rather 
reviving an old one) [18]. 

A. Gravitino Decay 

Since the gravitino has Planck suppressed couplings, 
it is never in thermal equilibrium in the early universe 
and depending on how much of it is produced at various 
reheating phases, it can come to dominate the energy 
density and ruin BBN just like the cosmological moduli 
can (in this paper we are assuming that the gravitino is 
not the LSP and that it decays). 

If this is the only reheating phase in the early universe, 
then the abundance is directly proportional to the reheat- 
ing temperature and solving the Boltzmann equations for 
the inflaton/radiation/gravitino system gives 

^3/2 - 2 X 10-12 __L (14) 

3/2 VlO^^GeVy ^ ' 

The BBN constraint on the temperature of reheat from 
inflation (T^^^) can be very stringent. For 7713/2 ~ 30 
TeV, y3/2 must be smaller than 2 x 10-^^ at 95% con- 
fidence level which implies T™-^ < 10^° GeV while for 



^ Note that in addition to the BBN constraints, the moduli is also 
constrained by WMAP measurements of isocurvature perturba- 
tions |23| . We leave it to future work to check the implications 
of these constraints on non-thermal production of dark matter 
from moduli. If dark matter is produced by the decay of a grav- 
itino, instead of a modulus, no isocurvature perturbations will 
be produced. 



smaller values (say TO3/2 ~ TeV) the bounds are more 
severe and T™^ needs to be as low as ~ 10^ Gev which 
is a very serious constraint on infiationary models [24 . 

Assuming that the temperature of reheat from infiation 
satisfies the constraints, the gravitino problem is revived 
in the presence of decaying moduli and a LTR [18' . In- 
deed assuming a branching fraction -B3/2 of the moduli 
(f) to '03/21 then from non-thermal production during the 
decay of 4>, we produce gravitini with 

where — can be obtained directly from the Boltz- 
mann equation for 0. Unless -B3/2 is tuned to be small, a 
low mass gravitino (again say less than 20 TeV) is ruled 
out in this scenario and the heavy moduli give rise to the 
gravitino problem. 

One way out is to assume that the moduli has sup- 
pressed couplings to the gravitino and that B3/2 is nat- 
urally small. On general grounds we expect this branch- 
ing ratio to be of order 0.01 — l[18] for » TO3/2. In 
[T9) . it was argued that it could be smaller due to he- 
licity suppression. In this case the decay width of the 
modulus to gravitino is suppressed from its total decay 

width T total ~ ^ to r^^^3/2 ^^1^. Given a hierarchy 
between and of order of Att'^ (as an example), 

we can get a branching ratio of order lO^-**^ ^ lO*-^*^^ 
which can be enough to evade the BBN constraints. Al- 
ternatively, if 7710 < 27773/2 then the branching ratio is 
drastically reduced due to phase space consideration. 

Given that the gravitino problem is so pervasive in 
many models of the early universe, it is tempting to as- 
sume that the abundance is not tuned (or diluted) to be 
small. The gravitino will dominate the energy density of 
the universe but if it decays and reheats prior to BBN, 
there will be no problem. To get this one needs a fairly 
heavy gravitino and one is lead to SUSY breaking pattern 
of the type of anomaly mediation or mirage mediation. 
In this case there is hierarchy between the gravitino mass 
and the LSP which we will parametrize 

7773/2 = Kuix . (16) 

In mirage mediation that we will discuss below, this hi- 
erarchy is of order k = Att'^. 

Now the gravitino itself will produce dark matter with 
a yield of 

V V STT^g* MpT3/2 (CTV) J 

where we have assumed a branching ratio of order 1 for 
gravitino decay into the LSP. The temperature of reheat 
for the gravitino T3/2 is determined in the same way as 



5 



for the cosmological moduli with a decay rate 

3 



(18) 



So T3/2 is given by Eq. Q replacing c ^ C3/2 and ^ 
TO3/2- In [H], it was shown that for MSSM dark matter 
candidates such as the Wino (with {av) ~ 10^^"'), the 
gravitino will overproduce dark matter unless its mass is 
higher than around ~ 10"^ TeV. A similar problem was 
pointed out in Moroi and Randall |2D] in the context of 
LSP production from modulus decay, where it was argued 
that a modulus mass of around 300 TeV gave ^lsp ~ 1 , 
while a lower value of f2 ^ 0.1 is obtained for even higher 
modulus mass. 

Assuming that ¥3/2 is large, the abundance of dark 
matter is given by the second factor in Eq. (17). The 
enhancement factor is given by 



6.14 X lO'^ 



3/2 



TeV 

mx 



1/2 



1 



(19) 



^3/2 



The branching ratio of the gravitino to LSP is essentially 
1 and C3/2 is constrained to be maximally around 1.5 
|25j (essentially, supersymmetry fixes the coupling of the 
gravitino to the supercurrent). For a hierarchy oi n — Att^ 
and a LSP mass around mx ~ 1.5 TeV an enhancement 
factor of Tf/Tr ^ 10* is obtained, which would require 
a cross-section of order {av) ~ cm? s^^ . As we will 
discussed more in the subsequent sections, it is possible 
to fit PAMELA with such a high cross-section although 
there is a definite tension with BBN constraints coming 
from dark matter annihilation. 

For a larger hierarchy, k = IGtt^, the enhancement 
factor is of order Tf /T^ ^ 10'^ and the cross-section must 
be {av) ^ lO^^^cm^s^^ which can give a very good fit 
to PAMELA as we will show below. 

Interestingly, the helicity suppression that can reduce 
^3/2 (and ^3/2) from the decay of the modulus is in 



general not enough to make the first factor in Eq. (17) 
smaller than the self-annihilation abundance. For the 
range of numbers used in this paper, we found that the 
branching ratio needs to be smaller than J53/2 ^ 10"'' for 
the first factor to be smaller while the helicity suppres- 
sion of [19 gives 10"'' — 10~^. This is interesting as this 
means that even when BBN constraints for the gravitino 
are evaded by having -B3/2 small enough, one still over- 
produces the LSP unless the cross-section is enhanced by 
the factor we calculated. Of course, if the modulus mass 
is small and the branching ratio is suppressed because 
of phase space consideration, then i?3/2 can be as small 
as we want. In that case the non-thermal production of 
dark matter will be dominated by the modulus and not 
the gravitino. Since they would have very similar tem- 
perature of reheat, we expect about the same enhance- 
ment factor although in principle the modulus could have 
a very small branching ratio to dark matter (unlike the 
gravitino) . 



To summarize, the non-thermal production of dark 
matter from the decay of the gravitino can give rise to 
the correct relic density if the cross-section is larger than 
the canonical value by a factor of 10^ — lO'*. This is in 
the high range of what is allowed by experiments but it 
could be the explanation for PAMELA as we will further 
discuss below. If the gravitino abundance is always very 
small, (even smaller than what is required to satisfy BBN 
constraints), then we can neglect its contribution to the 
dark matter relic density and instead look at cosmologi- 
cal moduli. By tuning the mass of this modulus one can 
get enhancement factor between 1 — 10**. In this case one 
must worry that the gravitino problem is not revived in 
this process by ensuring B3/2 is small enough. 



III. A CONCRETE MODEL 

In this section, we construct a successful model to im- 
plement the cosmological scenario oulined in the previous 
section. We study mirage mediation |^ to a U{1)b-l 
extension of the MSSM which appears often as a typical 
setting from the point of view of string phenomenology, 
and argue on general grounds that the dark matter is ei- 
ther the right handed sneutrino or the U{1) b-l gaugino. 
Further the dark matter is leptophilic since the dominant 
mode of annihilation is to the lightest of the new Higgs 
fields, and its subsequent decay mainly produces taus or 
muons by virtue of appropriate B — L charges. The cos- 
mological moduli and gravitino problems are addressed 
by the rather large cross section in such annihilation. 

In a construction of supersymmetry breaking vacua 
such as KKLT [27], the volume modulus T is stabilized by 
non-perturbative effects and then an uplifting mechanism 
(for example with anti-D3 branes) is used to obtain a su- 
persymmetry breaking vacuum. The exponential form of 
the non-perturbative potential leads to a small hierarchy 
between the moduli mass and the SUSY breaking scale 



Tot ~ TO3/2ln(Mp,/TO3/2) 



(20) 



with ln(Mp;/TO3/2) ~ 47r2. If the supersymmetry break- 
ing brane is sequestered from the visible sector, then T 
makes 0{Ft/T) ^ m^^/mT contributions to the soft 
terms and one sees that the modulus contribution to soft 
masses is comparable to the anomaly-mediated contribu- 
tion. There is thus a natural hierarchy of sparticle, grav- 
itino, and moduli masses given by 0{Att'^). Note that 
in general one might expect extra moduli with masses 
around the SUSY breaking scale to^ ~ ^3/2 and while 
they can be used and included in the discussion we mainly 
discuss the minimal mirage scenario with very heavy 
moduli in this paper. Setting the LSP at the TeV scale, 
we obtain the cosmology depicted in Fig. 1. 
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FIG. 1: The main thermal events in a scenario with a heavy 
moduh at 10"^ TeV, a gravitino mass at 10^ TeV and a LSP 
at the TeV scale. We assume that there are no suppressed 
branching ratios. Enough LSP (denoted X) is produced for 
self-annihilation to be important and the attractor solution 
for the abundance is reached at each phase transition. At each 
reheating phase, there is entropy production and the previous 
abundance of dark matter is diluted (by a factor roughly of 

( "r^id ) between the two phases 'new' and 'old'). The final 
answer in this particular set-up is to a good approximation 
simply given by the last decay Yx ~ m^t^T^t^- 



The MSSM with a B - L extension 



For gauge group SU{3) x SU{2)l x x U{1)b-l, 

C7(l)'g_^ triangle anomaly cancellation automatically 
implies the existence of three right-handed (RH) neutri- 
nos through which one can explain the neutrino masses 
and mixings |28j . Such extensions have been studied for 
a long time; this model has also been recently studied 
in the context of inflation [SHj, anomaly mediation [5U] . 
dark matter [TTl |3T], and leptogenesis The model 

contains a new gauge boson Z' , two new Higgs fields H'l 
and H2, and their supersymmetric partners. The B — L 
charge assignments arc shown in Table 1. The superpo- 
tential is 

W = W^MSSM + Wn + tJ.'H[H'2 + Wsojt (21) 

where Wjq is the superpotential containing RH neutrinos, 
and \j! is the new Higgs mixing parameter. Note that 



(22) 



where yD corresponds to Dirac Yukawa couplings. The 
new Higgs fields are neutral under MSSM charges, and do 
not have renormalizable couplings to quarks and leptons. 
The 1](X)b-l symmetry is broken by the VEV of these 



new Higgs bosons, v\ = {H[) and = {H2) 
Yukawas generate small neutrino masses. 



The Dirac 



Fields 


Q 




L 




N 


N" 


H[ 




Qb-l 


1/6 


-1/6 


-1/2 


1/2 


-1/2 


1/2 


3/2 


-3/2 



TABLE L The B - L charges of the fields. Here Q, L and 
TV represent quarks, leptons, and RH neutrinos respectively; 
while H[ and H2 are the two new Higgs fields. The MSSM 
Higgs fields have zero B — L charges and are not shown in the 
table. 



The new bosons have masses as follows. The U{\)b-l 
gauge boson Z' receives mass m?^, = {21 /4)g'^_^{v'i + 
V2), with gs-L being the B — L gauge coupling. There 
are three physical Higgs states: the lightest </> (not to be 
confused with cosmological moduli which we considered 
in previous sections) has mass to^ < m\, cos^ 2/3', which 
implies ni^ <^ mz' for tan /?' = v'2/v'i ~ 1. The other 
two Higgs states, <& and A, are heavy and have masses 
comparable to mz'- 



If our new Higgs have charges 1 and -1, we can al- 
low couplings: fijH'iN^N^. In this scenario, the Majo- 
rana masses for RH neutrinos develop through fij, once 
v[ 0. If we consider the case of B — L breaking at the 

TeV scale, then, the see-saw mechanism to generate small 

2 2 

neutrino masses requires ^ ^^^"^ '--^ 0.1 eV, which gives 
yfj ^ 10~^"/. We further note that / ~ 0.2 (assuming / 
to be the largest Majorana coupling) is large enough for 
radiative U{1)b-l breaking, and small enough to guar- 
antee that one-loop corrections to do not dominate 
over its tree-level bound. For this value of /, one obtains 
< 20 GeV which is needed to satisfy the anti-proton 
data from the PAMELA experiment. 

We can have either Z' or the right-handed sneutrino 
N as the LSP. The masses of the new Higgsinos are de- 
termined by fi' and they are assumed to be heavy. The 
lightest neutralino in the B — L sector lies predominantly 
in the Z' direction. 



The beta function coefficients for the B - 
given by 



L model are 



B-L 



61, 62, 63) = (51/4, 33/5,1,-3). 



(23) 



At the electroweak scale qb-l ^ 0.4, in order to achieve 
grand unification. Note that a normalization factor of 
•y/3/2 has been used for the B — L charges. The anoma- 
lous dimensions for the matter superfields (taking yB 
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much smaller than the MSSM Yukawas) are given by 



are 
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(24) 



If we introduce the Majorana couplings / in the model 
then 7jv<: and jh' will get — contributions in the above 
equation. We are not including / in our analysis just for 
simplification and even if we include these couplings, the 
overall conclusion remains unchanged. 



B. Mirage Mediation 



Mirage mediation is a mixture of modulus and anomaly 
mediation. In this scheme of mediation, the gaugino and 
scalar masses unify at an intermediate scale (called mi- 
rage scale) below the GUT scale. This scheme occurs 
quite naturally in warped compactification of string the- 
ory such as in [57] but we will not rely on any specific 
string theory construction. The soft parameters at the 
GUT scale are given by 



Ma = Mo-} 



^ijk 
TO? 



^3/2 

167r2 
_ TO3/ 2 

IGtt 

"^3/2 



bagl, 

2 (7* + 7j 



TO^^Mo.. 



Ik), 



(25) 



where Mg, Aij^., and TOj are pure modulus contributions, 
given as functions of the modulus T. Our conventions 



-3tr(r,2(Adj))+^tr(r,2(0,)) 

i 



li 



a jk 



! Ajjk 

Mo ' 



(26) 



where the quadratic Casimir C^{(pi) = (iV^ - l)/2Ar 
for a fundamental representation ipi of the gauge group 
SU{N), Cf (0,) = qf for the U{1) charge of cjii, and 
'^ki VikiUjki is assumed to be diagonal. To set input pa- 
rameters, we define the ratios 



MoHMpi/my^)'' 



Mo'' 



Mr 



. (27) 

where a represents the anomaly to modulus mediation 
ratio, while and Ci parameterize the pattern of the 
pure modulus mediated soft masses. 

The input parameters in RG running are 



Mo, a, Oi, Ci, tan/3. 



(28) 



where one could also choose TO3/2 as an input in place of 
a. 

In terms of brane constructions in type IIB, if the 
matter fields live on the entire worldvolume of the Dl 
from which visible sector gauge fields originate, then 
fli = Ci = 1 while if the matter fields live on intersections 
of D7s, then Oi = q = 1/2,0 [55]. Compactifications 
with dilaton-modulus mixing, realized, for example, in 
type IIB by the presence of gauge flux on the Dl , can 
easily lead to other positive values of a, c, and a. In typ- 
ical compactifications, TO3/2 is set by appropriate choice 
of flux contributions to the superpotential 



Ad 



aT 



(29) 



with ln{A/Wfiu^) ^ in^ and m3/2 W^flux- 

For typical values of the model parameters, the gaugino 

mass runnings are shown in Figure 2. 

The ratio of the gaugino masses to ^, : TOBino : TOwino : 

TOgiuino is obtained as 



(1 : 1.2 : 1.8 : 3.6). 



(30) 



Evidently, the Bino is heavier 
mediation to the MSSM, the hj 
Bino for a < 1, with the Hig 
ing with increasing a. This is 
OijCi, and for tan/3 = 10 [25] . 
in the B — L extension, since 
MSSM Yukawas get negative 
slightly raising their low enerj 



than Z' . In usual mirage 
jhtest neutralino is mostly 
;gsino component increas- 
true for various values of 
This conclusion changes 
the beta functions of the 
contributions from qb-l, 
ly values. This difference 
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FIG. 2: Running of gaugino masses in the B — L model. 
The thick gray, sohd, dotted, and dashed-dotted lines are 
the gluino, Wino, Bino, and Z' respectively. We have used 
ai=Ci=l, tan/3 = 10, mg/a = 77TeV, Mq = 2.5 TeV (a = 1) 



feeds positively into the beta function of the MSSM fi pa- 
rameter, thus lowering its low energy value compared to 
pure MSSM. In principle, this would mean that the Hig- 
gsino component of the lightest neutralino in the MSSM 
sector would begin to dominate for slightly lower values 
of a. 

In Figure 3, we plot the RG evolution of the sfermions. 
For m3/2 = 77 TeV, Mq = 2.5 TeV, and q = = 1 (a = 
1) one sees that the scalars are heavier than the B — L 
gaugino. The RH sneutrino is lighter than the MSSM 
sfermions, due to the fact that in the case of sleptons 
we have contributions from MSSM gauge couplings in 
addition to the U{1)b-l gauge couplings. 

We can make the right handed sneutrino even lighter 
by choosing Ci, appropriately. We show one such example 
in Figure 3 (the solid line at the bottom) and in fact, the 
right handed sneutrino can be the LSP of this model. For 
this case, we take, 7713/2 ~ 200 TeV, Mq = 5 TeV, Ci, — 

— 0.3, Ci = Oi = 1 for all other particles (a — 1.3), one 
obtains the right-handed sneutrino as the LSP with mass 
around 1.5 TeV while the gauginos and other scalars are 
around ~ 3 TeV or so. The cross section in this case is 
around ^ 10"^'^ cto^s"^. 

So we conclude that cither the B — L gaugino or the 
right-handed sneutrino can be the LSP in this model. 



C. Explanation of the observed anomalies in 
cosmic rays 

In order to explain the recent cosmic ray data, we need 
electron-positrons in the final states of LSP annihilation. 
In this model, the LSP {Z' or the sneutrino) annihilates 
into light Higgs bosons (from the B — L sector) which 
then decays into a pair of taus predominantly jTT] . 

The taus then decay into electron-positron pair. A re- 
cent analysis of the data showed that [11^ in order to 
explain the excess by using rs, we need an enhancement 



factor of 10'^ for the annihilation cross-section. The an- 
nihilation cross-section does not have any p-wave sup- 
pression. The typical value of the LSP mass that fits 
the data in this model with this enhancement factor is 
about 1-2 TeV. The LSPs annihilate to hghtest Higgs (0) 
of the B-L sector, whose mass is controlled by the VEVs 
of the new Higgs fields. For comparable VEVs, i.e. for 
tan/3' « 1, it can be very small without any tuning of 
the soft masses in the Higgs sector. We can choose this 
mass to be between 0(1) GeV and 20 GeV in order to 
be in complete agreement with the anti-proton data. For 
2mT- < TO0 the dominant decay mode is to t~t+. If 
we assume the ip mass to be >20 GeV, then the Br of 
(/)—!■ 66 is about 1/7 of ^ rr due to the B-L charges. 
In this case the anti-proton data is still satisfied up to 
a factor of 2. However, the computation of anti-proton 
flux involves a large theoretical uncertainty |33] . In our 
model both small and large values of cj) are allowed, and 
as we have already discussed, a small Higgs mass requires 
smaller values of tan /?'. If to^ is slightly less than 2mr, 4> 
can decay either to cc or fJ-^ fJ-^ with comparable branch- 
ing ratios. It is possible to reduce the mass further to 
be below 2toci and make final state the dominant 

decay mode. 

We now discuss the two options for the LSP in this 
model. 

1. If Z' gaugino is the LSP, then from the cosmol- 
ogy analysis we find a typical enhancement ^10*^ 
gives the correct relic density. This value of the en- 
hancement follows from Eq. ( [T9| ) and the fact that 
the gravitino mass and the LSP mass are typically 
correlated by a hierarchy n ^ \n.{Mpi/ra^/2) ^ 47r^. 
One can lower the enhancement by raising Mg, but 
in our case we cannot raise Mq too much since we 
want the gaugino mass to be 1.5 — 2 TeV. 

An enhancement factor of lO'* can have problem 
with BBN [31]'', where it is claimed that the en- 
hancement factor should be less than 10^"'^. If this 
enhancement factor is somehow accommodated by 
the BBN data, then it is possible to generate such 
a large cross-section by having an annihilation fun- 
nel oi B — L gauginos into a pair of the 0, the 
lightest boson in the B — L Higgs sector via the s- 
channel exchange of the cj), <& (heavy Higgs). The 
S-channel resonance of this process enhances the 
cross-section to the required value. While a 1.5 
TeV LSP has been fitted with PAMELA data for 
an enhancement 10^ , we note that astrophysical un- 
certainties (for example, a choice of isothermal DM 
density profile instead of NFW in the halo func- 
tion can easily give rise to a factor 2-5 uncertainty. 



* The most stringent bounds coming from ^He/D give the upper 
bound on the cross section for non-thermal dark matter annihi- 
lating to tau to be 



' 3 X 10 ■^•^ erne's 
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and a factor of 2 uncertainty in the energy loss co- 
efficient of positrons [35 ) could provide a fit with 
enhancement 10*. 

The best option for the LSP is the new sneu- 
trino (-/V). In this case, the hierarchy between 
the gauginos and the gravitino remains ~ 47r^, 
but for small values of in the mirage media- 
tion input parameters, the sneutrino can be made 
much fighter than the gauginos. Thus, the hierar- 
chy between the sneutrino LSP and the gravitino 
becomes k ~ IGtt^ and a lower cosmological en- 
hancement of 10^ (hence annihilation cross-section 
around 10~^^cto^s~^) is obtained from Eq. (19 1. 
An enhancement factor of lO'^ allows us to fit the 
PAMELA data. The BBN bound of [34^ is also 
satisfied without any difficulty. On the particle 
physics side, this cross-section can be obtained with 
and without the heavy Higgs annihilation funnel. 
The annihilation amplitude is proportional to the 
gauge boson mass which appears in the N*N(/) ver- 



tex ^. The dominant channel is N*N 



via the 



s-channel exchange of the <&, the t, u-channel 
exchange N, and the contact term |7Vp(/)^. The 
s-channel Z' exchange is subdominant because of 
the large Z' mass (as required by the experimental 
bound on mz'). The sneutrino annihilation into 
vP final states is at least an order of magnitude be- 
low the (f>(j) final states. Other fermion final states, 
through s-channel Z' exchange, have even smaller 
branching ratios (these fermion-anti-fermion final 
states are p-wave suppressed). 

Since the cosmological enhancement is sufficient to 
explain the PAMELA data, we do not need any en- 
hancement due to Sommerfeld effect. Sommerfeld 
enhancement requires [11] the model parameters to 
be tuned very accurately, and this can be easily pre- 
vented. 

One interesting aspect of sneutrino LSP is that they 
can be probed in direct detection experiments [11]. The 
direct detection cross-section for sneutrino-nucleon scat- 
tering is mediated by Z' exchange and the cross-section 
can be quite large. To explain the PAMELA data in our 
model we need to use large Z' gauge boson mass (if we 
do not use the heavy Higgs annihilation funnel) and con- 
sequently the direct detection cross-section is reduced. It 
is interesting to note that if instead we use thermal dark 
matter with Sommerfeld enhancement in our model, a 
smaller Z' mass is needed to explain the dark matter 
content (again, if we do not use the Higgs annihilation 
funnel). Thus the direct detection cross-section in the 



This vertex arises from Vd D 
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FIG. 3: Running of right sneutrino (A*'), left selectron and 
right selectron masses in the B — L model. From top, the 
solid, dotted, and dashed-dotted lines are the sneutrino, left 
selectron and right selectron masses, respectively. We have 
used a,=Ci = l, tan/3 = 10, mg/a = 77TeV, Mo = 2.5 TeV 
(a = 1). The solid line which does not go through the mirage 
point corresponds to N for Cy=0.3. 



case of thermal dark matter is much larger, 
this feature in Figure 4. 



We show 



Therefore, combining the direct detection result with 
PAMELA results it is possible to distinguish the cosmo- 
logical enhancement from the Sommerfeld enhancement. 
If we choose the annihilation funnel to satisfy the dark 
matter content, we can allow smaller values of Z' and the 
direct detection cross-section can be larger. 
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FIG. 4: Sneutrino-nucleon scattering cross-section as a func- 
tion of sneutrino mass. The top line corresponds to the ther- 
mal case where we need Sommerfeld enhancement to explain 
the PAMELA data. The bottom line corresponds to the non- 
thermal case as discussed in this work. The correct relic den- 
sity is satisfied for both cases. 
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IV. DISCUSSION 

In this paper, we have studied non-thermal dark mat- 
ter, the moduH and gravitino problem in the light of the 
PAMELA data. As a phenomenological model, we have 
considered a U{1)b-l extension of the MSSM where su- 
persymmetry breaking is mediated by mirage mediation. 

We have found that the final decaying particle that 
non-thermally produces LSP may be either a cosmologi- 
cal modulus or the gravitino. Cosmological moduli typ- 
ically produce gravitino, and the decay of either at a 
temperature above BBN tends to overproduce dark mat- 
ter. A larger annihilation cross section for dark matter 
can naturally ease this overproduction problem. 

We have shown that it is possible to solve the mod- 
uli/gravitino problem in the B — L model with mirage 
mediation. The natural hierarchy between LSP and grav- 
itino/moduli in mirage mediation allows the gravitino to 
decay above BBN, while maintaining an LSP in the 1 — 2 
TeV range. Moreover, a large enhancement of the anni- 
hilation cross section of lO^"** is readily obtained in this 
model, which solves the overproduction problem and can 
fit the PAMELA data. 

The LSP can be either the new B ~ L gaugino or the 
right-handed sneutrino. Both of these annihilate to the 
light Higgs of the new sector. This Higgs primarily de- 



cay into tau for to^ > m(2r) due to the B-L charges. 
For the B — L gaugino LSP 1.5 TeV, an enhancement 
~ 10^ is calculated cosmologically from the decay of the 
gravitino. This enhancement does not do well with BBN 
constraints. For the sneutrino LSP ^ 1.5 TeV, a cos- 
mological enhancement of 10'^ is calculated from a larger 
hierarchy between LSP and gravitino. This can be ob- 
tained by an appropriate choice of mirage mediation pa- 
rameters. This enhancement explains the recently ob- 
served anomalies in cosmic rays and is allowed by the 
BBN constraint. The sneutrino LSP has interesting con- 
sequences for direct detection experiments. In fact, in 
this case it is possible to distinguish between models of 
non-thermal dark matter and models with thermal dark 
matter that utilize Sommerfeld enhancement in certain 
regions of parameter space. 
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